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Abstract 



The effects of competing quadrupolar- and spin-glass orderings are in- 
vestigated on a spin-1 Ising model with infinite-range random p-spin inter- 
actions. The model is studied through the replica approach and a phase 
diagram is obtained in the limit p — > 00. The phase diagram, obtained 
within replica-symmetry breaking, exhibits a very unusual feature in mag- 
netic models: three first-order transition lines meeting at a commom triple 
point, where all phases of the model coexist. 
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The mean-field theory of Ising spin glasses is quite well understood at the present ||T]J 
Since the pioneering solution of Parisi for the infinite-range interaction Ising spin glass, 
the so-called Sherrington-Kirkpatrick (SK) model 0, a wide variety of spin-glass systems 
has been investigated within the replica approach In particular, this method was 

used to solve a generalization of the SK model with the inclusion of p-spin interactions 
1^, and it was found that the p oo limit leads to the solution of the random-energy 
model, introduced earlier and solved by quite different methods 0. Since then, many 
other p-spin interaction models have been studied, motivated by the fact that they are 
tractable within mean-field theory, for arbitrary values of p, thus rendering it possible to 
analyse both p — > 2 and p — > oo limits 0. Another fundamental aspect of such models is 
their striking connection to real structural glasses 0. Also, it is usually feasible to study 
the dynamical properties of those models, making it possible to gain some insights in the 
important ageing phenomena presented by random systems [§,0. 

An important class of random systems, with many physical realizations, is that of 
orientational glasses [TU,rT|. These systems are usually described in terms of an assembly 
of discrete spin variables with quadrupolar random interactions. Recently, a p-spin inter- 
action orientational glass model was investigated by the replica method |T2|, and it was 
shown that Parisi's replica-symmetry-breaking (RSB) scheme could be applied success- 
fully: indeed, it was found that, in the limit p ^ oo, the low-temperature behavior of such 
model may be properly described through a single-step RSB approach. This conclusion 
was shown to be true for another class of models which describe multipolar glasses ||13 |. 

Despite all the above-mentioned efforts, much less has been studied on magnetic mod- 
els where different kinds of disorder are present. This represents a very common situation 
in physical systems, opening a wide variety of new problems to be investigated. For 
spin-1 Ising variables, a simple model including pairwise dipolar and quadrupolar random 
interactions led to interesting behavior, with a competition between quadrupolar- and 
spin-glass orderings |]T1 . 



In the present letter we investigate a spin-1 Ising model including both dipolar and 
quadrupolar random p-spin interactions. We consider an infinite-ranged interacting sys- 
tem, consisting of spins described through the Hamiltonian, 
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'H — ^ Jij^i^-ipSij^Si^ ■ ■ ■ Sip ^ Kij^i^...ip(Si-^Si^ ■ ■ ■ SipY , (1) 

l<ii<---<ip<N l<ii<---<ip<N 

where each spin variable can assume the values 0, ±1. Both couplings, Jji...jp and Ki^...ip 
are quenched, independent and identically distributed random Gaussian variables, with 
zero means and variances J'^p\/{2NP^^) and K'^p\/{2NP^^), respectively. 

It should be mentioned that a spin-1 Ising spin-glass model with p-spin interactions. 



under a single- ion anisotropy field D was already studied by Mottishaw [T^; such a 
model is identical to the present one only for D = K = 0. Futhermore, the resulting 
phase diagrams of these two models share some common features, e.g., they both present 
three distincts phases. However, whereas Mottishaw's model exhibits one first-order and 
two continuous critical frontiers, we show that the phase diagram of the present model 
displays no continuous lines, being characterized by three first-order critical frontiers, 
which meet at a triple point. 

Applying the replica method for the model defined through Eq. (1) and following 
standard procedures, we get the free-energy density 

- (3f = lim -Gn{qab, Aafe, Qab, lab, Ra, Ca), (2) 

where 

Gniqab, Xab, Qab, lab, Ra, Q = ^ Y.iJ' <fab + K^Qlb) + ^^^^^ ^ E ^1 

^ a^b ^ a 

- J E iXabQab + labQab) " E ^^Ra (3) 
^ a^b a 



-{- In Tr exp 



\ Y.iXabS^'S' + labiS'^S'f) + E Ra{S^) 
^ a^b a 



with P = {kBT)~^, and a,b = 1 . . . n denoting rephca indices. The quantities {Xab, lab, ^a) 
represent Lagrange multipliers, allowing to fix the set of order parameters {qab, Qab, Ra)- 
By demanding G„ to be stationary with respect to each of those parameters, we get the 
equilibrium conditions 

qab = {S^S'), Qab = {{S'^S')'), Ra = {{S'^)'), 

(4) 

^ab — 2 ^ab , lab — 2 ^ ab , Sa ~ 4 ^^a ■ 
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Throughout most of this letter, we will be interested in the limit p ^ oo. This is 
justified by the fact that p- interaction models usually exhibit the same qualitative behavior 
for finite values of p {p > 2). Besides that, finite values of p require a substantial amount 
of numerical work, whereas in the limit p — > oo most of the calculations may be carried 
analytically, based on the fact that the energies of distinct configurations are uncorrelated 

As a preliminary approach to the problem, let us consider the replica-symmetric (RS) 
solution, i.e., qat = q, Qab = Q, Ra = R, Kb = A, 7ab = 7, = in terms of which 
the equilibrium equations become 

with 

Q = {^l)xy , Q = {^l)xy , R = {^2)xy ■ (6) 

In the equations above, {{...))xy stand for / ^|^exp(-^-|^) (...), and 

2e^ siBh{y/Xx) 2e'' cosh(v^x) 
^1 = ^ ? V2 = ^ , (7) 



where 



S = -^^ + ^ + ^y^ Z = 1 + 2e'' cosh(v^: 



x] 



One may easily see that the order parameters Q and R never vanish and that Eqs. (|5|) 
and (|^) always present a trivial solution with g = A = 0, for arbitrary values of p; herein 
we identify such a solution with a quadrupolar-glass (QG) phase . 



In the limit p ^ oo the QG solution represents the only acceptable solution associated 
with the parameter q, since any solution with g 7^ is unstable, similarly to what happens 
in the corresponding spin-| model 0; however, there are two simple solutions for the 
parameters Q and R, as we describe below. The first one is given by 7 = ^ = 0, and 
Q = = 4/9, in which case the free-energy density becomes 

f = -kBT\n3, (9) 



4 



corresponding to an entropy per spin s = fcs ln3; we shall refer to the phase described by 
such a solution as a quadrupolar-glass 1 (QGl). A second solution can easily be found 
with Q = R = 1, 2'y = jS'^K'^p, and 4^ = f3'^{J'^ + K'^)p, in such a way as to yield a 
free-energy density, 

f = --f--ksTln2. (10) 

This solution, which we call quadrupolar-glass 2 (QG2), presents an entropy per spin 
which becomes negative for ksT/J < kBTc/ J = ^-^^ = 0.6005 ■■ ■. For T > Tc, both 
solutions are stable, the former one (QGl) presenting a lower free energy at high tem- 
peratures. As the temperature is lowered, we find a first-order transition line, where the 
free energies of those solutions coincide; this line is independent of K and is given by 
knT^lJ = — , ^ = 0.7852 •■ -. It is important to mention that we have also found 
other solutions, all of them being completely unstable. Therefore, in the limit p oo the 
RS solution leads to the phase diagram exhibited in Fig. 1, with two quadrupolar-glass 
phases, QGl {q = 0,Q = 4/9, R = 2/3) and QG2 {q = 0,Q = R = 1), separated by a 
first-order transition line. 

Since the RS Ansatz leads to a QG2 solution which becomes unstable at low tem- 
peratures, one must carry on with a RSB procedure. In analogy with the spin-| cor- 
responding problem [^,0 and other p-interaction orientational glasses |]12|,|13|, one may 



see that it is sufficient to consider a single-step Parisi RSB scheme for the present prob- 
lem. This is achieved by grouping the n replicas into n/m blocks of m replicas each. 
Order parameters with replica indices a, b in the same block take on certain values 
{(lab = Qi, Qab = Q15 ^ab = -^1, = 7i), distinguished from those with replica indices in 
different blocks {qab = Qo, Qab = Qo, ^ab = Aq, 7a6 = 7o), whereas the single-replica-index 
parameters are considered in the RS approximation {Ra = R, ^a = 0- usual [^, in 
the limit n ^ the parameter m becomes a continuous variable in the interval [0, 1]. For 
an arbitrary value of p, the free-energy density becomes 

= — [{m - l)g^ - mq^] — [{m - l)g^ - mQ^o] Rp 

+ i [(m - l)(Aigi + 7iQi) - m(Aogo + loQo)] + 



1 f°° dxodxi 



+ -/ / =^exp(-^^)lnZ(a:o,Xi), (11) 
m J-oo J-oo 27r ^ 



where 

zM = ^rr 'AM-ii^m.-r. (12) 

m J-oo J-oo ZTT Z 

and 

B = (/^xo + ^Ai - \oyo)S + (Vt^xi + VtT^I/i + C - y - y )^'- (13) 

Prom the above free-energy density functional we can obtain several solutions, including 
those already described within the RS approximation. From now on, we will restrict 
ourselves to the limit p ^ oo. Our analysis indicates that there is only one new solution 
within a one-step RSB that is physically acceptable, i.e., 50, Qq < l,qi — Qi — R — 1. 
This solution presents a a free-energy density independent of T, given by 



/ = -^{P + K^)\n2 , (14) 

with a zero entropy. It corresponds to the low-temperature phase, where all spin variables 
are frozen completely at random, each of them in one of the states Si — ±1, like in a 
spin- 1 Ising spin glass. Such a solution will be associated with the quadrupolar-spin-glass 

(QSG) phase. 

Since we have obtained three phases (QGl, QG2 and QSG) and their respective free 
energies, we can draw the phase diagram of the model within the RSB approach (see 
Fig. 2). We adopt the standard thermodynamic criteria, i.e., whenever two or more 
solutions are stable, the correct phase is defined as the one with the lowest free energy. 
Besides the first-order transition line separating phases QGl and QG2 (which remains 
a line independent of K at ksTi/J = 0.7852- ■ ■), we find two new first-order transition 
lines: one represents the coexistence of phases QG2 and QSG and is given by 

whereas the other one corresponds to a coexistence of phases QGl and QSG, 



The three phases coexist at a triple point, whose coordinates are given by 
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kBTt/J 



1 



0.7852 ■ ■ ■ , KtlJ 



21n2 - ln3 



0.2713- •• , 



(17) 



2^111(3/2) 



2y^ln21n(3/2) 



where all three lines merge together, in accordance with the famous Gibbs phase rule ||16|| . 

In conclusion, we have solved a disordered spin-1 Ising model with p-spin interactions 
through the replica method. Considering the limit p — oo, we have verified that the 
Parisi Ansatz is suitable to determine correctly a phase diagram with genuine first-order 
transitions, each of them accompanied by a latent heat and exhibiting discontinuities on 
the respective order parameters. Each phase is described by the solution which is a global 
minimum of the free energy. A very uncommon feature in magnetic models has been 
detected, i.e., a triple point where all three phases coexist. An important aspect to be 
explored is the dynamics of this model, mainly along the first-order transition lines; this 
issue is currently been investigated. 
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FIGURES 



FIG. 1. The phase diagram within the RS approximation. The quadrupolar- glass phases, 
QGl and QG2 (defined in the text), are separated by a first- order transition line. Throughout 
the gray region the solution QG2 becomes unstable, presenting a negative entropy. 

FIG. 2. The phase diagram within Parisi's RSB procedure. All phase boundaries are 
first-order transition lines. At the triple point (black dot) the three phases coexist. 
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1 - de Araujo et al. 




Fig. 2 - de Araujo et al. 



